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We analyse the stability of global 0(3) monopoles in the infinite cut-off (or scalar mass) limit. We 
obtain the perturbation equations and prove that the spherically symmetric solution is classically 
stable (or neutrally stable) to axially symmetric, square integrable or power-law decay perturbations. 
Moreover we show that, in spite of the existence of a conserved topological charge, the energy barrier 
between the monopole and the vacuum is finite even in the limit where the cut-off is taken to infinity. 
This feature is specific of global monopoles and independent of the details of the scalar potential. 



I. INTRODUCTION 



o 
o 
o 

Global monopoles have been investigated for years as possible seeds for structure formation in the Universe |PL||. 
Although they appear to be ruled out by the latest cosmological data 0], their appearance in condensed matter 
-and other- systems and their peculiar properties make them worthy of investigation. These objects have divergent 
energy, due to the slow fall-off of angular gradients in the fields, which has to be cut-off at a certain distance R (in 
practice, the distance to the nearest monopole or antimonopole) and has two important consequences, in particular 
for cosmology. First, the evolution of a network of global monopoles is very different from that of gauged monopoles, 
as long-range interactions enhance annihilation to the extent of eliminating the overabundance problem altogether 
0. Second, their gravitational properties include a deficit solid angle Q|, which makes them rather exotic. 

The stability of global 0(3) monopoles has been the subject of some debate in the literature J5H7]]. In this paper 
we try to settle the issue by: a) analysing the axial perturbation equations in the limit where the cut-off is taken to 
infinity, and b) proving that the energy barrier between the monopole and the vacuum (meaning the extra energy 
required by the monopole to reach an unstable configuration that decays to the vacuum) is finite. It is somewhat 
surprising for different topological sectors to be separated by finite energy barriers, but in this case it is a consequence 
of the scale invariance of gradient energy on two dimensional surfaces (r = constant), and therefore independent of 
' the details of the scalar potential. 

II. THE MODEL 



in 
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We consider the simplest model that gives rise to global monopoles, the 0(3) model with lagrangian: 

1„ 1 



L = ^$ Q cW - ^A(|$| 2 - r? 2 ) 2 a = 1,2, 3. (1) 



$ a is a scalar triplet, |$| = V ( E >a( I )a and /a = 0, 1, 2, 3. The 0(3) symmetry is spontaneously broken to 0(2), leading 
to two Goldstone bosons and one scalar excitation with mass m s = \/2Xrj. The set of ground states is the two- 
sphere |$| = rj and, since 7T2(5 2 ) = Z, there are field configurations with non-trivial topological charge. One such 
configuration with unit winding is the spherically symmetric monopole, 

® a =Vf(r)- , (2) 
r 

where /(0) = and f(r — > oo) = 1. Its asymptotic behaviour is f(r — > 0) ~ or, a ~ 0.5 and f(r —* oo) ~ 
1 — 1/r 2 — 3/2r 4 , as can be seen from the e.o.m. of f(r) 

/" + - f (/ 2 - !) = ■ ( 3 ) 

The two parameters (77, A) appearing in the lagrangian can be absorbed by the rescaling <I> a — > $ a = < f >a /?7, x^ — » 
x^ = \J Ar/ 2 x M , which amounts to choosing 77 as the unit of energy and the inverse scalar mass as the unit of length 
(up to a numerical factor). Note, however, that the energy of a configuration with non-trivial winding such as (|^) 
is (linearly) divergent with radius, due to the slow fall-off of angular gradients and has to be cut off at r = R, say. 
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Unlike r\ and A, the (rescaled) cut-off is an important parameter which could affect the dynamics of solutions with 
non-trivial topology. Dropping tildes: 

E= td^ a d^ a + ~(|$| 2 - l) 2 . (4) 
Jo 2 4 

Since the energy diverges, Derrick's theorem does not apply in this case, and in JtJ it was shown that the global 
monopole is stable towards radial rescalings. On the other hand the question of stability with respect to angular 
perturbations has led to some discussion in the literature after Goldhaber Q pointed out that the ansatz 

"J 1 = F(r,0) sm0(r,e) COS If 
$ 2 = F(r,9) sin 6(r,9) simp 

$ 3 = F(r,6) cos S(r, 9) , (5) 

which describes axially symmetric deformations of the spherical monopole (0) , leads to the following expression for 
the energy after a change of variables y — lntan(6>/2): 

E = C d< ^ I-oc d V Jo" dr \ [Pi + r 2 sech 2 (y)p 2 ] 
Pl =F 2 +F 2 [sm 2 (0) + (9 y ) 2 } 
p 2 =F^ + F^9 r ) 2 + \{F 2 -l) 2 (6) 

with F r = d r F, etc. The term in brackets in pi is identical to the energy of the sine-Gordon soliton, so translational 
invariance in y implies that configurations with 

F(r, y) = /(r), tan (0(r, y)/2) = e a+ «, £ = const (7) 

have the same energy as (Q) (which corresponds to £ = 0). On a given r =const. shell, the effect of taking £ — > oo is 
to concentrate the angular gradients in an arbitrarily small region around the north pole. When the gradient energy 
is inside a region of size comparable to the inverse scalar mass, it is energetically favourable to "undo the knot" by 
reducing the modulus of the scalar field to zero and climbing over the top of the mexican hat potential. Unwinding 
is estimated to occur at a critical value of £ (say £q) whose dependence with r far from the core is logarithmic, 
£o w lnr + const. Numerical simulation on individual shells closer to the core gives £o ~ «o + &olnr with slowly 
varying bo ss 1 and oq ~ —1.3. Unwinding is expected if £(r) > Co( 7 ")- 



III. ESTIMATION OF THE ENERGY BARRIER 

As explained in 0,0, the shift (R) creates a tension pulling the monopole core, and the apparent unwinding (which 
starts in the inner shells) is only a manifestation of the core's translation. In order to stop the motion of the core, 
we consider a hybrid configuration such that the monopole core remain unperturbed and the unwinding occur in 
the outer shells. This is achieved by taking £ = for, say, r < r\ 1 a string- like configuration for r > and some 
continuous interpolation in between. One such configuration (see figure Ha) would be: 

!0 r < n 

c(l-a) n<r<r 2 (8) 
a + b ln(r) r 2 < r 

with a > ao, & > &0j an d c given by continuity of £(r). 

If the energy (i.e. mass) inside r 2 is large enough there will be no appreciable motion of the monopole core (because 
the tension of the string is constant ~ Air). For instance, simulations with ri = 3, r 2 = 6, 6 = 1 show different 
behaviour depending on c: for c < 2.8 the core translates, but for c > 2.8, the unwinding happens far from the core, 
which remains fixed. In the latter case, the decay of the string can also be understood as a monopole-antimonopole 
pair creation with the new monopole appearing at r = R and the antimonopole appearing near r — r 2 . The equations 
solved in the simulations were 

□ $ a + $ a (|$| 2 _l) +7 $ a = , (9) 
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a) 




FIG. L a) Plot of the scalar field <E> a in the configuration given by eq. ([]) in the text, for r\ — 3, ri = 6, a = 0.9, 6 = 1. 
The configuration is axially symmetric, b) The result of numerical integration using a) as initial condition. Potential energy is 



shown in greyscale at different times. After annihilation of the monopole-antimonopole pair near z — 0, 
for some time until all the energy is radiated away (not shown; note the longer time elapsed between t = 



the system fluctuates 
:2.5 and t = 4.5). 



with a dissipative term added to make the integration faster. Simulations using different 7 show no appreciable 
difference. The equations were integrated using cylindrical coordinates (p, <p, z) in a 200 2 grid using explicit Runge- 
Kutta method with step size control ||. The output can be seen in fig. [ljj, where the potential energy is plotted at 
various times, confirming that the configuration (|J) is unstable and decays to the vacuum. 

We will now show that the extra energy required to reach this unstable configuration (||) from the spherical monopole 
configuration (||) is finite. 

Consider the set of configurations with F — f(r) and £ = £(r). From (||) the difference between the energy of any 
such configuration (-£[£]) and the energy of the spherically symmetric unperturbed monopole (i?[0]) is: 

y/ 2 M£ r V)/K(r)] 

/+00 
dysech 2 (y)sech 2 (y + e(r)) . (10) 
-00 

/[£] is bell-shaped: from a maximum value I[0] — 4/3, it rapidly falls to zero for |£ | > 5 as ~ 16(|£| — l)exp(— 2|£|). 
r* is the radius at which £(7%) = The first integral is clearly finite. The second can be estimated using the 
asymptotic form of and is negligible for £ = a + lnr even in the limit R — > 00, where one gets ~ 16 7re _a ^e _l ** . 

Moreover, there is a continuous path in configuration space connecting the configurations (^|) and (||) such that 
E[£] — E[0] remains finite along the whole path: first increase £ in the outer shells r > r* using Goldhaber's deformation 
(Cr = 0) until £ reaches £*; then adjust the radial dependence to match (||). 

Thus, we have shown that the extra energy required by the monopole (Q) to go over the energy barrier and 
decay to the vacuum is finite even as R — > 00; moreover, since the monopole energy grows with R, the ratio 
(E[i] - E[0])/E[0] -> as R -> 00. 



E[$ - E[0] = 
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FIG. 2. The functions: a) f(r)/r, f r (r); b)Vi(r), V 2 (r). 



IV. STABILITY TO SMALL PERTURBATIONS WITH AXIAL SYMMETRY 

We now turn to the classical stability of (^|) by considering small perturbations parametrised by 

F = f(r)+S(t,r,y) 
tan(fl/2) = (l + C(t,r,y))e«' 

Neglecting quadratic terms gives sin 8 as sin 8 + £ sin 8 cos 6, cos 8 s» cos 6* — £ sin 2 6*. Introducing X = /£, 

$ : « (/ sin + 5 sin + X sin cos 0) cos ip 
$ 2 ps (/ sin 8 + S sin 8 + X sin cos 8) sin <£> 
$ 3 w (/ cos 6» + 5 cos 6> - X sin 2 0) 



(11) 
(12) 



(13) 



which shows that the correct boundary conditions are that sm8(5 + Xcos8) should vanish on the z axis. Note 
that 5(0) and X(0) need not vanish. An infinitesimal translation of the monopole in the z direction corresponds to 
6 = —f r (r) cos 9, X — f(r)/r, and both f r (r) and f(r)/r tend to w 0.5 as r — > (see fig ^a ). There is no zero mode 
associated with global rotations since these have been factored out in the ansatz (|J). 
As usual, the perturbation equations 



= r 2 U5 + 25 + r 2 [3f 2 - 1]6 + 2X y - 4Atanhy 
= scch 2 y r 2 UX + sech 2 y[r 2 {f 2 - 1) + 2]X + 2tanhy X a - 25 y 



(14) 



reduce to an eigenvalue problem in lu 2 for perturbations of the form <5 = e lult 5(r, y), X = e llAjt X(r,y). Eigenfunctions 
5, X with negative eigenvalues oj 2 < correspond to instabilities. Dropping hats and defining u = tanhy 



RiS + d u [{u 2 - l)d u 5] - 2d u [{u 2 - l)X] = uj 2 r 2 S 
R 2 X + d 2 [{u 2 - 1)X] - 28 U S = LU 2 r 2 X 



(15) 



where R%, i?2 are radial operators (see fig |^b): 



Ri = ~d r (r 2 d r ) + V t (r), 
Vx(r) = r 2 (3/ 2 (r)-l) + 2 
V 3 (r) = r^r) - 1) 



i = 1,2 



(16) 



Using Legendre polynomials and changing variables: x '■— <9 U [(1 — u2 )X] = J2xi( r )Rii u )i $ — J2^i( r )Pi( u ) the 
equations for different values of I decouple, giving 
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i?iA ; + x 2 Ai + 2x X , 
R2X1 + x 2 xi + 2xA, = lo 2 t 2 xi 



(17) 
(18) 



where we introduced x — ylQ + 1) and A; = x5i. In order to get ( fL8[ ) we multiplied the X eqn. in ( |l5| ) by (1 — u 2 ) and 
differentiated w.r.t. u, so there may be spureous solutions; in particular, I = corresponds to angular perturbations 
that are singular on the z axis. These are not physical, and will be discarded. But if there is no solution of (17 1^) 
with negative w 2 there will be no instability in the original problem ([l4|). 

Our task is to find the solution to (^,^) with the minimum value of to 2 over all admissible perturbations and 
all I > 1. We know one 1 = 1 solution, the translational zero mode (Ai = y/2f r {r), xi = -2/(r)/r). Goldhaber's 
deformation (xi = f(f), Ai = 0) is also I = 1 but is not a solution of (|l7j , |l8| ), and it can be shown that there are no 
instabilities with Xi( r ~ * 00 ) ~ const. 

Let us first consider normalizable perturbations. Note that, for each I, the equations (|l7],|l8|) can be obtained by 
functional variation from 



E t = J dr [r 2 (A 2 . + X 2 ) + (Vi + x 2 )A 2 + (V 2 + x 2 ) X 2 + ±xA X ] = lo 2 J 



drr 2 \A 2 



The lowest value of to 2 can be found minimising Ei over normalized functions, J drr 2 [A 2 



X 



(19) 



1 and over all 



I > 1. However the minimum must be in the I = 1 sector since, for all i > 1 and for given A,%, (Ei — E\) is a sum of 
squares with positive coefficients: 



Ei — Ei 



dr[A+(A + X ) 2 +A-(A-x) 2 ] 



(20) 



where Ai t ± 



2 ± 2{x - V2)]/2 > Ofor I > 1. 



In order to investigate the Z = 1 sector, we rewrite £1 using arbitrary functions G(r), i?(r) and K(r): 
Bh. = J °° dr [(rA r + f A) 2 + (r Xr + fx)' + 2^2 (A + X x ) 2 + (^1 + 2 - + G r - ^) X 2 



T2 



2V2K 2 



H r -K)A 



G(r)A 2 (r 



1 00 
lo 



H{r) X \r) 



(21) 



Choosing G(r) = — r 2 9 r lnAi, H(r) = — r 2 <9 r lnxi, K(r) 2 = — Ai/xi, the coefficients of A 2 and x 2 in ( ^l| ) vanish 
identically by virtue of (||). For large r, G and ~ r, thus, i?i > for all functions (Ai, xi) that decay faster than 
as r — * 00. This proves that all normalizable perturbations have w 2 > 0. The above argument and a host of 
numerical simulations strongly suggest that non-normalizable perturbations are at best compatible with to 2 = 0, as 
can be verified directly from the equations for perturbations that fall to zero like a power of r, but in this case we 
have no analytic proof. 



V. DISCUSSION 



In this paper we have derived and analysed the axial perturbation equations of 0(3) monopoles and proved that, 
contrary to statements in the literature, 0(3) monopoles are perturbatively stable (or neutrally stable) to infinitesimal, 
axially symmetric, normalizable (or power-law decay) perturbations. We have also proved that the energy barrier 
between topological sectors is finite, irrespective of the details of the scalar potential. This feature is specific of global 
monopoles; global vortices in two dimensions, whose energy grows logarithmically with radius, and gauge monopoles, 
whose energy is finite, are separated from the vacuum by an energy barrier growing linearly with the cut-off. But the 
energy of global monopoles is dominated by two-dimensional (angular) gradients far from the core, and these can be 
deformed with no energy cost. 

One would naively expect thermal fluctuations with KT ~ E[£] — E[0] to cause the monopole to decay. As far as 
we know, this effect is not seen in "cosmological" numerical simulations of global monopole networks, but perhaps 
this is not so surprising. First, we worked in flat space. Second, the range of scales introduced by the expansion of the 
Universe forces drastic approximations on cosmological simulations; in particular, the sigma model approximation, 
which is widely used, sets the field on the vacuum manifold everywhere, so unwinding and pair creation events are not 
resolved by the grid. Finally, cosmological simulations do not include thermal effects, and it has been shown that full 
thermal simulation across the phase transition [ pl)[ gives qualitatively different results. Our results provide further 
evidence that a more careful analysis of global monopole networks may be required. 
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